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The eigenvector expansion developed in the preceding paper for a system of damped linear oscil- 
lators is extended to critical points, where eigenvectors merge and the time-evolution operator TL 
assumes a Jordan-block structure. The representation of the bilinear map is obtained in this basis. 
Perturbations eATL around an M-th order critical point generically lead to eigenvalue shifts ~ e 1 ^ 1 
dependent on only one matrix element, with the M eigenvalues splitting in equiangular directions 
in the complex plane. Small denominators near criticality are shown to cancel. 
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I. INTRODUCTION 

The preceding paper [l| [hereafter referred to as I and 
equations therein as e.g. (1-2.3)] developed an eigenvector 
expansion for a broad class of systems with N ohmically 
damped coupled oscillators, also applicable to interact- 
ing quantum systems. The key concept is a bilinear map 
(i/>,</>), under which the time-evolution operator TL is 
symmetric: (tp,TL<f>) = (Hi/), 0), thus allowing concepts 
familiar from conservative systems to be transcribed. 

In I, it is assumed that (a) the eigenvectors fj of TL 
are complete and (b) {fj,fj) ^ for all j. These are 
violated (together) only at critical points (with measure 
zero in parameter space) , where eigenvectors merge. This 
case is rather more technical, and TL takes on a Jordan- 
block (JB) structure. Some results are known in a con- 
tinuum limit 0, but this paper gives a self-contained 
account, putting the concept of critical damping into a 
general framework, in which projections etc. can still be 
performed using the bilinear map. 

Since merged eigenvectors no longer span the whole 
space, a basis has to be constructed — the well-known Jor- 
dan normal basis :3| . Here, we need to represent the bilin- 
ear map (viz., the metric g..) as well and verify that the 
resulting TL.. is symmetric (Section[nJ|. Then time evolu- 
tion (Section |ll F|) and perturbation theory (Section |III|) 
are developed in this basis. The latter also permits a 
discussion of small denominators near criticality. 

An N = 1 example (Section III A in I) already illus- 
trates criticality, and also explains its name. Section ITVl 
presents higher-order examples, including perturbations 
around them. We end with a discussion in Section^ An 
extensive account can be found in Ref. 0. 

We stress that when eigenvalues merge (degeneracy), 
either the eigenvectors merge as well (criticality), or they 
could remain distinct (level crossing) — though the latter 
is exceptional in a sense to be made precise in Section llVl 



II. JORDAN BLOCKS 



For any finite matrix TL. & Jordan normal basis is 
known to exist in general In a revival of interest 
in non-hermitian and non-diagonalizable systems [f| , JBs 



increasingly have the attention of the physics community. 
We briefly recapitulate some standard results in linear al- 
gebra, and address the remaining fruitful object of study: 
the consequences of JBs for TL of the specific form (1.2.3), 
corresponding to the well-motivated dynamics (1. 1.1) and 
leading to the bilinear map (1.2.9). 



A. Basis vectors 

Suppose there are v independent eigenvectors fj with 
eigenvalues uij. If v < 2N, the dimensionality of phase 
space, the fj can be augmented into the Jordan normal 
basis {fj^iKjKvflKnKMj, obeying 



(TL — Ulj) fj, n — fj, n -l , 



(2.1) 



with fj t —i = so that fjfi = fj. Separating posi- 
tions and momenta, fj >n = {fj, n ,fj, n ) T - The span of 

{fj, n } n=o 1 f° r nxe d j is called the JB at Uj, with size Mj\ 
in it, TL has the Jordan normal form 



TL'. = 
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(2.2) 



with respect to the normal basis. The block structure 
{(ujj,Mj)} is completely specified by TL, but vectors from 
blocks with identical eigenvalues may be mixed, as al- 
ready shown by degeneracies in the diagonalizable case. 
Until Section III Dl however, all Uj are assumed distinct. 



(2.3) 



Then, the only arbitrariness consists of 



fj 



(*k fj,n — k 



with Co ^ 0, leaving H2.2|l invariant and generalizing fj i— > 
cfj for simple ujj . This freedom will be exploited below. 

When JBs occur, J(uo) = det(W. — u>) has a multiple 
zero. But this is not sufficient, since it could also cor- 
respond to level crossing without JBs. This is why the 
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formalism is developed without reference to J — in con- 
trast to the case of continuum models (Section VI in I). 



B. Bilinear map 

Although the Jordan normal form itself is standard, 
we need to consider the bilinear map (fj.n, fj> ,-«')■ Using 

(EH) on (Tift ,71; Jj',7%' ) = (fj,n,^ifj',n>) y ields 

fj',n') + \ fj,n-l,fj',n') 

= u 3'(fj,n> fj',n>) + (fj,n>fj',n'-l) j ( 2 -4) 



It is now straightforward to verify that 7i.. is 



for n < Mj — 1 and n' < My 



1. 



Orthogonality for j ^ j' is proved by induction with 
respect to n + n' . The case n = n' = is the stan- 
dard one as in (1.2.11). On each side of l|2. 4J1 . the sec- 
ond term vanishes by the induction hypothesis, leaving 



(ujj- 



')(fj,n, fj',ri) — 0, and completing the induction. 



Next consider j = j'; with M = Mj and (n, n') = 
(fj,n,fj,n>), <E3I> gives (n-l,n r ) = (n,n'-l), leading to 



(n, n') = A n+Ti 



(2.5) 



For n < M—2, we have A n = (—1,71+1) = [again as- 
sociating incomplete eigenvectors with (fj,fj) = 0]. If 
Am-i would also vanish, then /^o would be orthogo- 
nal to every basis vector, which is impossible [cf. below 
(1.2.16)]. We can choose Am— i = 1 and A n = for 
M < n < 2M—2. To show this, first perform a transform 
(|2.3|) with Cfc = C()Sko, under which Am-i * CqAm-i = 1 
for some cq ^ 0. Further transforms = Sko + c n 5kn 
for n = 1, ...,M— 1 take Am+u-i = [n, M— 1) i— > 
Am+u-i + 2c„ = by a choice of c„, while Ajvr+m-i 
with m < n (taken care of in previous steps) are not 
affected (neither, of course, arc A n = for n < M—l). 
All these orthogonality relations are captured by 

(fj,n,fj',n') — §jj'8 n +n',Mj-l , (2-6) 

and completeness can be written as 

<t> = ^2 fj,n {fj,M t - l-n, <t>) ■ (2-7) 



Equation (|2.6I) can also be stated as a representation 
for g.. in (one block of) the Jordan normal basis: 
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(2.8) 



n.. 






uji 



to 




0/ 



(2.9) 



a symmetric Jordan-type matrix. 

Therefore, the duals defined by [cf. (1.2.19)] 

/'" nj\,, >,../, u 



ensure 



(f j,n \fj',n') =Fj' S n n < 

^ = E^« (/ i,n i->. 



(2.10) 

(2.11) 
(2.12) 



the latter in an obvious shorthand. In this case, T><p can- 
not be calculated without first resolving it in terms of the 
basis; contrast the non-critical case in I. Of course, the 
(unique) dual to any basis {e/j} can always be calculated 
by inverting the 27V x 27V Gramm matrix (efc|ej) Q with- 
out any reference to the bilinear map, but this matrix is 
not sparse for {e^} = {fj, n }- The above has bypassed the 
general inversion in that all manipulations below l|2.5|) 
take place within a (typically small) block only, made 
possible by 1)2. 4|) . Cf. the remarks in I, Section II E. 

Together with our normalization, l|2.1|l implies that the 
counterpart to (1.2.14) for conjugate blocks reads 

/_ J -,» = ±t^(-) n / J >. (2-13) 

where the overall sign should be the same within each 
block. For a block with imaginary frequency (a so-called 
zero-mode, see I) which does not cross with another 
block, we can set j = 0, and (|2.13|l becomes a symme- 
try of the basis vectors. The treatment of level crossing 
follows the one in I and will not be repeated. 

Incidentally, in most cases only one nontrivial JB 
(say j) is formed at a given parameter value. Then, the 
fj >n can be constructed without having to solve the lin- 
ear system H2.ll . First, determine fj\ now if> = [g-fj]* is 
guaranteed to have (ip, fj) ^ 0. Orthogonalize ip with re- 
spect to all the other eigenvectors [which does not change 
(ip,fj) 7^ 0]; the result is f^M -\- After the fj iTl with 
lower n are obtained from 12.1|l '. the normal basis at ujj 
can be biorthogonalized as above. 



C. Example 

Consider a single oscillator as in I, Section III A, but 
at the critical k — fc* = 7 2 . The eigenvector is /.o = 
c (l,- 7 ) T . Then/,! =c H/7,0) T + Cl (l,-7) T . 

One verifies = and A\ — Cq, independent of C\\ 
we take Co = 1. Finally A2 = 2(ci — i/7) = if c\ = i/j. 
With these choices, f t o = (1, — 7 ) T and f .i = (0, — i) T . 
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D. Level crossings 

When levels cross in conservative systems, orthogonal- 
ity is only a matter of choice. Here the situation is more 
subtle still, since blocks of arbitrary size may cross. Sup- 
pose that the normal basis has uj\ — ui^ = . . . = wl ISl- 
The corresponding vectors [their number given by the 
order of the zero in J (to)] form an array, e.g., 

/l,3 /2,3 

fl,2 h,2 f3.2 (2 M) 

Jl,l J2,l J3,l 

fl,0 /2,0 /3,0 fifi 

each of the L columns corresponds to a JB, obeying l|2.1[) . 
While degeneracy is doubly exceptional (as will be seen in 
examples) , we give the rigorous treatment for the record. 
We need to ensure 

(fj,n,fj'n') = Sjj'Sn+n^Mj-l ■ (2-15) 
It is immediately shown, in the usual way, that 

(/i,«,/i',n0=4+n'> ( 2 ' 16 ) 

reducing the number of conditions to be checked. To 
proceed, we use the freedom 

L n 

hn -EE ^'ff.n-k , (2.17) 
j'—l h—rn 

where m = max(0, Mj — Mji). The mixing between 
blocks is the new feature here, as is apparent from the 
case of two degenerate trivial blocks. The lower limit 
k = m signifies that large blocks can always mix into 
smaller ones, but not vice versa. For instance, in l|2.14[l . 
the only freedom in /^o and /2,o consists in mixing these 
two vectors among each other. 

We start by finding a normalizable "top" vector. If 
M% < Mi, then fi Mi—i can be normalized by the argu- 
ment below 1)2.50 . If there is no single largest block, con- 
sider the quadratic form (</>, (H—cu j) Ml_1 4>) m the span 
of all vectors [ m the example l|2.14[) . j = 1,2]. 

This form cannot vanish identically, for else the corre- 
sponding eigenvectors /^o would be orthogonal to the 
entire space by l|2.16[l . Take one vector on which the form 
is nonzero, relabel it as /i Mi— 1> and redefine the other 
h" necessary so that they still span the same sub- 
space. Subsequently, construct the associated lower fj n 
by iterating "H — uij. Clearly, the basis transformation of 
this paragraph is of the form (|2.17l) . 

In the block j = 1, one now has precisely the situation 
of Section lll Bl consequently, this block can be biorthog- 
onalized as described there. It remains to orthogonalizc 
the other blocks with respect to this first one. Then rela- 
tion i|2.15|) will hold for j = 1 and any j' , which will not 
change when the blocks j > 2 are mixed among them- 
selves upon repeating the whole procedure. 



For orthogonalizing the blocks j > 2 with respect to 
the j = 1 block, completing one iteration of the con- 
struction, it suffices to consider the top vectors fj,M -i- 
when their associated block is reconstructed by iterating 
TL — ojj, all vectors in it will be orthogonal to {/i,„} by 
virtue of l|2.16|l . By the same token, (fj^Ms-x, fi,n) = 
already if n < Mi — Mj. For n > Mi — Mj, this bilin- 
ear map can be made to vanish by /^m^-i — /i,M,-i + 
^ L j+ „_ Ml /i,Af 1 -i-n- Thus, the needed mixings are pre- 
cisely those which conserve the block structure and are 
allowed by (|2.17() : e.g., in l|2.14|l . in the course of the 
above we do not mix /13 into fz.2- 

Iterating the procedure until the number of blocks L is 
exhausted, we finally achieve l|2.15|l in general. Clearly, 
all of the above goes through as well if level crossing 
occurs at several frequencies. 

E. Sum rules 

Separating (|2.7|) into coordinates and momenta leads 
to four sum rules, which can be written in terms of the 
coordinates alone. Using the shorthand n' = Mj — 1 — n, 

° = 2/j,n«>/j,n' (2.18a) 

i=J2 + ( 2 - 18b ) 

= E [( W jfa n + 2w j'/j,n-l + fj,n-2)®fj,n' 

+ i{uj 3 f 3 . n + 4»-i)®(T/ iin /)] (2.18c) 
= X)/j,»®(l7 J ,„')- (2-18d) 

3." 

Here, a ® b is the matrix with elements a(a)b((3). Wc 
have verified these in, e.g., the examples of Section llVI 

F. Time Evolution 

Since TL is not diagonal, the basis vectors' time depen- 
dence is slightly complicated. The defining idtfj, n (t) = 
UjfjA 1 ) + fj,n-i{t) [fj,n(t=0) = fj,n] is solved by 

n 

fjA*)=Y l C l (u j ,t)f J , n -i , (2.19) 
1=0 

= t^t e -^ = if^e"*-*]^ . (2.20) 

The evolution of a general initial state 0(t=O) = <fi fol- 
lows by simply putting fj iTl 1— > fj, n (t) in l|2.7|l . More 
formally, the retarded Green's function thus reads 

Mj-l 

m = mJ2 E hn(t)(fn-) ■ (2.21) 

j n=0 
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Since 6(t)Ci(ujj, t) > i(u> — Uj) under Fourier trans- 
form, this leads to 



Mi -i 



W^EE^CT(FI') - (2-22) 

j n=0 1=0 y 31 

readily verified to solve (H — uj)Q(uo) — —il. 



III. JORDAN-BLOCK PERTURBATION 
THEORY 

A. Lowest order 

Although perturbation theory has been given in I, the 
situation at a critical point is different, exhibiting inter- 
esting features not found in conservative systems. For 
simplicity consider only the case without level crossing. 
Under a perturbation of a critical Ho given by l|2.2ll in 
the Jordan normal basis, 



H = H + eAH 



a JB of size Mj to leading order generically behaves 
as follows, (a) The eigenvalue splits into Mj different 
ones, shifting in equiangular directions in the complex 
frequency-plane and at the same rate. The directions for 
e > bisect those for e < 0. (b) The frequency shifts go 
not as e but as e 1 ' M ' . (c) The shifts depend on only one 
element of the Mj x Mj matrix AH within this block. 

These features are already seen in the trivial N = 1, 
Mj = 2 example of I, Section III A. For k = fc* + e, 
ujj = — 17 ± y'e, approaching the critical point along the 
real direction for e > (slightly underdamped) and along 
the imaginary direction for e < (slightly overdamped). 

To derive the above properties, consider the eigenvalue 
equation. For the moment, focus on one block; inter- 
block couplings will be added in Section IIII Bl without 
difficulty. From (|2.2I) . for its lowest-order determinant, 
H — lu is effectively represented by 



/{jJj—LO 1 

H'.-lj = 



\ < o 

where in the determinant 



\ 




(3.2) 



u)j~u)/ 



= i-fiMj-i 



\AHf j ,o) = (fj, ,AHf j , ) 



(3.3) 



multiplies the l's, and all other elements of AH are neg- 
ligible because they multiply another small quantity of 
at least order uj 3 — uj = — Alu. Consequently, in 

J(lu) = (-l) M J [J H + eJ 1 H + e 2 J 2 H + ---] , (3.4) 



one has Jq(lu) — (uj — uij) j and Ji{u)j) — — £. Setting 
J(uS) = (-Aw) M > + (-lp*-^ = 0, one finds 



Auj k = X(k 







_ 2iirk/Mj 



(3.5) 
(3.6) 



A is any fixed choice of the root measuring the magnitude 
of the shift, displays the equiangular behavior, and 
k = 0, . . . , Mj— 1. A change in sgne costs a phase e i7r / Mj , 
causing the directions to bisect the original ones. 

Incidentally, for the conservative case with degenera- 
cies, ll.'i.2t would not have the l's and the crucial term £ 
does not appear at this low order. 

The eigenvectors can be expressed as 

M } - 1 

fk= ]T f 3 , n T n k ; (3.7) 

n=0 

upon setting uj = ujp. in l|3.2|l . one solution is 



T n k = (AC fe ) r 



( 3 - 1 ) [so (T-^n = (ACfe) - ™^], with normalizatic 



(f k ,f k )=Mj(\t k ) 



M,--l 



(3.8) 



(3.9) 



The above applies to any H cast into Jordan normal 
form. But in the present case, H is given by (1.2.3); 
in particular, we assume that AH does not affect the 
coupling to the bath, so Ar = as in (1.4.1). Then, 



(3.10) 



in terms of the coordinates only. In the example of I, Sec- 
tion III A, if we let K = k* + e, i.e., AK = 1, then (|3~TU|) 
gives £ = 1 so l|3.5(l gives Auj k = \^e Z7rk , as expected. 



B. Higher orders 

To deal with higher-order corrections, we rewrite 

H = H + eAH = H' + eAH' , (3.11) 

shifting the term e£ to the unperturbed part, so that H3.7(l 
diagonalizes H' . In other words, (fi^-^AH'fjfi) = 0. 

The matrix elements of AH' in the split basis l|3.7|l 
follow by transforming those in the Jordan normal basis, 

[AH'f y = {T- V ) k n (r< n \AH'fj, n ,)T n ' k , 

Mi-l 
3 n,n'— 

(3.12) 

Because the term (n'—O, n=Mj — l) has been removed 
to H' , the lowest power of A is \ 2 ~ Mj . Together with 
e = X Mj , the leading correction is 0(A 2 ), one power of A 
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higher than the effect due to £ treated in Section Till Al 
it can be handled by the standard perturbation theory 
as in I. Each higher-order correction will involve an extra 
matrix element oc A 2 , divided by u>k — ^k' °^ \ gi vm g one 
overall power of A per order. Inter-block matrix elements 
of AH' can be handled in the split basis as well. 

In fact, there is not much to be gained by further ana- 
lytic treatment; one can simply diagonalize the relevant 
Mj x Mj block of H numerically and deal with inter- 
block interactions in the resultant basis, using the non- 
degenerate perturbation theory of I. 



C. Non-generic perturbations 

The perturbation of a JB is said to be non-generic 
when £ = (AK)jj = 0. The simplest way to proceed is 
to expand (|3.4I) around u>j, 

= (-l) M U(u>) 
= (Au) M i + e[Ji(^) + AujJ[(ujj) + ■■■} 
+ e 2 [J 2 (uj j ) + •••] + ••• . (3.13) 

In this case, J\(ujj) = 0; assuming there is no higher- 
order nongenericity, i.e., J[(u)j) ^ 0, (|3.13|) gives 

= Aw[(Aw) M ^ 1 + eJ' x {ujj) + •••] + ••• . (3.14) 

Thus we see that (a) one state is unshifted to lowest or- 
der, and (b) the other states split like a generic JB of 
order Mj — 1, viz., with |Aa>| cx e x ^ M '~ l ") and splitting 
in Mj — 1 equiangular directions. [If Mj = 2, the e 2 J2(uij) 
term is of the same order and must be retained.] This 
situation will be seen in some of the examples in Sec- 
tion \W\ The exercise analogous to IpT ^ -lpO )) and (|3~TU|) 
shows that J[(u>j) = — 2£', with 

£' = (/,-,!, AHf jt0 ) = fj,i (a)AK(a, 0)f jfi (fl) ■ (3.15) 

Not much is gained by formally pursuing nongeneric 
perturbations any further, since in applications it is again 
preferable to simply diagonalize the (typically small) JB. 



and also (£5) as e~ luJkt = e^'* Ei>oH A Cfc07 z! = 
J2i>o(Kk) l Ci(uj j ,t), we get 

l>0 3 n,n'=0 k 

(3.17) 

where m = n-\-n'-\-l + l — Mj can be negative. But 
by (|3.6|) . J^k van ishes unless m = (mod Mj). Since 
the above is only valid to leading order anyway, we only 
keep m = terms, so Ylk'" ~ an< ^ l|3.17(l reduces 
to 4>{t) = Y, n /j,n(*)(/ J> l0) by Thus, the small 

denominators cancel (i.e., negative powers of A do not ap- 
pear): while the contributions per mode to 4>(t) are large, 
their sum remains finite. Moreover, the latter agrees with 
(|2.21(l , obtained in the Jordan normal basis at u>j [to low- 
est order; in general, <f>{t) smoothly depends on system 
parameters as these go through their critical values]. 

Similar arguments prove the cancellation of small de- 
nominators in other physical quantities, or for other split- 
ting patterns ( Section IIII C(l . These exactly parallel the 
discussion of excess noise (over the standard Schawlow- 
Townes value) in lossy laser cavities In fact, only 

the noise per mode is enhanced while the sum over all 
(non-orthogonal) modes has no excess contribution, con- 
sistent with the fluctuation-dissipation theorem. Only 
perturbation theory requires special care: while a near- 
critical mode cluster shifts other eigenvalues by a finite 
amount, the cluster's modes themselves are highly sensi- 
tive to perturbations ( Section IIII Bfl — indeed, no cancel- 
lation is expected in the properties of individual modes. 



IV. EXAMPLES OF HIGHER-ORDER 
CRITIC ALITY 

We construct some examples with JBs of size M > 2. 
It suffices to consider an N = 2 system, with 



An 


k 12 \ 


, r = 2( 711 


7l2\ 


Ul2 


&22 / 


V7i2 


722/ 



involving 6 parameters in general. 



D. Small denominators near criticality 

The above resolves the "small-denominator problem" . 
Time evolution near a critical point u>j is given by 

£ \fk,fk) 

where the prime denotes restriction to one cluster of 
modes fk, which are close to merging at uij as param- 
eters are tuned. Already in part I it has been argued 
that then (fk, fk) —> while (fk, <fr) is finite, apparently 
causing 4>(t) to diverge. However, by using (|3.7|l - (|3.9|l 



A. Fourth-order JB 

To find an M = 4 JB, let J(uj) have a 4th-order zero: 

det(H'. -w) = (uj + iaf (4.2) 

(necessary but not sufficient, since equality of eigenvalues 
could also indicate level crossing; the further conditions 
will be considered later). The position of the root has to 
be on the negative imaginary axis (or else there would 
be another root at —to*). Without loss of generality, we 
henceforth set a = 1, and the a-dependence can eventu- 
ally be restored by scaling K i— > a?K , V i— > aT, u> i— > auj. 
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Eliminating the momenta in (Ti' . — aj)<p = 0, i.e., set- 
ting uj = — i in (1.2.7), gives 



(K -T + !)(/) = 



(4.3) 



Generically, there is only one independent satisfying 
l[4.3|) . hence only one eigenvector <fi, hence an M = 4 JB. 

Equation (|4.2|l leads to four conditions for the coeffi- 
cients of u> 3 , . . . , uj°, which after some simplification give 



knk 



22 



k 



12 



1 

7n + 722 = 2 

2 



fell + ^22 + 4(7n722 ~ 7l2) = 6 
&11722 + &227U - 2^12712 = 2 



(4.4) 



The general solution for K involves two parameters: 

_ fe y cosh x sinha; \ , , 

I sinhx e~ y cosh.x) ' \ ■ ) 

solving the first of the four constraints. For any x, y, 
the three remaining equations determine the three 7«. 
Positivity of K is guaranteed, while T > requires 
cosh x cosh y < 3. A simple choice (which satisfies this as 
an equality) is sinha; = — 2, e 2y — 5, giving 



K = 



r = 



4 




(4.6) 



The Jordan normal basis is constructed to be 

/,o = >/2i (1,1,-1, -1) T 
/,i = ±V2(-1,1,3,1) T 
f, 2 = |V2*(-1,-1,5,-3) T 
/, 3 =^V2(-1,1,-1,-3) T , (4.7) 



normalized as in (|2.6|) . These vectors are arbitrary up 
to one overall sign, and exemplify the formalism in Sec- 
tion [H] In particular, the alternation of real and imag- 
inary basis vectors is prescribed by (|2.13|) . here realized 
with the lower sign. The duals follow from H2.10(l as 

/' =iV2i(5,3,l,-l) T 



f 1 = W2 (-1,3,1,1)' 



f' 2 
f< 3 



|V2i (1, -1, 1, -1) T 
v / 2(-3,l,-l,-l) T , 



(4.8) 



and are readily verified to obey l|2.11|l . 

It is possible to have more than one 4> satisfy (|4.3|l . For 
a 2 x 2 system, this requires K — T + I = 0, which to- 
gether with (|4.5|l fixes 7^ . When put back into the three 
remaining equations in l|4.4|l . these lead to one additional 
condition cosh x cosh y — 1 , which has only the solution 
x = y = — the trivial case of two independent but iden- 
tical oscillators, each generating an M = 2 JB at the 
critical point. Apparently, with N = 2, one cannot have 
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FIG. 1: Eigenvalues for an M = 4 JB split by fcn 



with e = n eg, n = 0, 1, 2, 



(a) 60 = KT 4 ; (b) e = -10 



Crosses (circles) denote numerical (perturbative) values. The 
nearly equal spacing shows that the shifts are oc e 1//4 . 



a crossing between M = 3 and M = 1 blocks. Thus, level 
crossing occurs only exceptionally, under the additional 
condition that a minor determinant of TC . — uj vanishes. 
The JBs of the two oscillators can be mixed, illustrating 
the subtleties encountered in Section lllDI However, tak- 
ing complex superpositions will in general break the sym- 
metry (|2.13() . generalizing the discussion below (1.2.14). 

It is instructive to consider perturbations around the 
M = 4 JB gUJl-lO. First, let ku 1-* fen + e, and 
evaluate the lowest-order result (|3.5|) (with £ = —2) for 
various e; these are shown by the circles in Fig. ^ illus- 
trating the properties below (|3.1|l . The exact numerical 
eigenvalues, shown by the crosses in Fig. ^ demonstrate 
the accuracy of perturbation theory. Second, consider 



the more general ki 



efitj, with fiij = O(l). By 



expanding J(oj) for the model (|4.1|l as in (|3.4I) . we find 



•M - *) = ( fc 22+l-2722)Mn + (&ii+1-27ii)M22 
+ 2(27i 2 -fei 2 )Mi2 • 



(4.9) 
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FIG. 2: Eigenvalues for an M = 4 JB split by the non- 
generic perturbation kij ^ kij + eflij, /til = 1, /*i2 = — §, 
M22 = 2; e = n 3 e , n = 0,1,2,.... (a) e = HT 4 ; (b) e = 
— 10 -4 . Crosses (circles) denote numerical (perturbative) val- 
ues. Three of the eigenvalues split at 120°; the nearly equal 
spacing shows that the shifts are oc e 1 ^ 3 . The fourth, nearly 
unchanged, eigenvalue is shown in greater detail in the inset. 



root has been scaled to —i. As before, one obtains 

fcn&22 — kf 2 = b 

711+722 = (3 + 6)/2 
kn + k 22 + 4(7n722-7i 2 2 ) = 3(1 + b) 
fcn722 + fc 22 7ii - 2fci27i2 = (1 + 36)/2 . (4.10) 

We work in the eigenbasis of T (cf. Section III B and 
note [3], both in I) from the outset, and can solve for 
the remaining parameters in terms of 711 and b. This 
involves rationals only for, e.g., 6 = 4 and K = f (Is 1 !)' 
r = (01), where for variation we took a V > 0. 
The basis vectors are found to be 

/i,o = e OT/4 #(2,-4,-2,4) T 

fi.i = e~ l7r/4 # (-19, -22, 43, -26) T 

/1.2 = eW4 M (~ 221 ' ~ 78 ' 525 ' 43 °) T 

/ 2 ,o = e"/ 4 ^(8,-l,-32,4) T . 



(4.11) 



It is seen that the j = 1 (j = 2) block obeys (|2.13(l with 
the lower (upper) sign; apparently, not much can be said 
about this sign in general. The dual vectors are 

/i.o = e "/ 4 ^|| (801, -352, 221, 78) T 



,-»7r/4 yT5 



2.1) 



„mt/4 yTS 
15 

a i7r/4 yTB" 
45 



^(-71,-48,-19,-22)^ 
^(-10,0,-2,4) T 



(-16,-3,-8, 1) 1 



(4.12) 



again up to an overall sign. 

Next change ku i— > fcn+e, giving £ = j^i; Fig.|3]shows 
the eigenvalues emanating from the JB (conventions as 



before). For fcj 



efiij, the choice /Jn 



-2, 



/ii2 = ^, ^22 = 1 makes £ = 0, eliminating the leading 
term, while £' = 1 in the next order. Figure 0] shows the 
eigenvalues; to leading order, one is not shifted, while the 
other two split like a generically perturbed M = 2 block. 



For the parameters of (|4.6p . this vanishes if /in = 1, 
/ii2 = — I , /*22 = 2, yielding a non-generic perturbation 
[£ = in 13-lOj l]. Figure |3 shows the eigenvalues vs e 
(crosses and circles as in Fig. To leading order one 
is not shifted, while the others split like a generically 
perturbed M = 3 block (i.e., ~ e 1 / 3 , smaller than the 
typical shift ~ e 1 / 4 ) — in quantitative agreement with the 
perturbative (|3.14(l and Ij3.15|l . where presently £' = i. 



B. Third-order JB 

Next consider a third-order JB, mainly to show that 
odd M 1 is allowed. For N = 2, this is achieved by 
setting J(ijj) — (uj + i) 3 (uu + ib), with 6^1. The triple 



C. Two second-order JBs 

So far, the examples only involve JBs with imaginary 
frequencies. This is not necessary: consider two second- 
order JBs at lo — — i ± b, by setting J{lo) = (oj + i — b) 2 x 
(uj + i + b) 2 . As before, this leads to 

k U k 22 - k 2 12 = (1 + b 2 ) 2 

7n + 722 = 2 
fcn + k 22 + 4(711722 - 7i 2 2 ) = 6 + 2b 2 
fcn722 + fe27n - 2/ci27i2 = 2(1 + b 2 ) . (4.13) 

Guided by Section IIV Al without loss of generality we 
proceed in the eigenbasis of T. With 712 = 0, l|4.13|l can 
be readily solved for the remaining parameters in terms of 
b and 711. Again we consider the marginal case 711 = 2, 
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FIG. 3: Eigenvalues for an M — 3 JB split by fen fen + e, 
with e = n 3 e , n = 0, 1, 2, . . .. (a) e = 10" 4 ; (b) e = -1(T 4 . 
Crosses (circles) denote numerical (perturbative) values. The 
nearly equal spacing shows that the shifts are oc e 1//3 . 



and a final simplification occurs for b = |, in which case 

all kij turn out rational: K = f ( _gj ~2° ) and r = ( o o ) ■ 
The basis vectors read 

f\,o = if (3-6i, — 3— 6z, -ll+2i, -5+10i) T 

/n = ^ (15+30i, -15+30*, -23-74i, 7+14i) T , 

(4.14) 

the conjugates f-\ n following from i|2.13[l . Their duals, 
/ 1)0 = ^§ (-46-37i, -14-7i, -30-15i, -30+15i) T 



f 1 ' 1 = *f (22-i,-10+5i,6-3i,6+3i)' 1 



(4.15) 



obey H2.11|l not only with the vectors in l)4.14[l . but also 
with their conjugates f-i, n . 

Figure El shows the eigenvalues when ku i— > fen + e, 
where in this case £ = — (9+12i)/32. (Because of sym- 
metry, only the right-half w-plane is shown.) 

Finally, consider (|4.13|) for variable b, keeping T = 
diag(4,0) fixed. One readily solves for K(b), which 
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FIG. 4: Eigenvalues for an M = 3 JB split by fey kij+e/j,ij, 
with /in = —2, /i22 = 1, Mi2 = 1/2 chosen so that the 
perturbation is non-generic, and e = n 2 eo, n = 0,1,2,.... 
(a) eo = 10 -4 ; (b) eo = — 10 -4 . Crosses (circles) denote nu- 
merical (perturbative) values. Two of the eigenvalues split 
at 180°; the nearly equal spacing shows that the shifts are 
oc e 1//2 . The third eigenvalue is nearly unchanged, and is 
shown in greater detail in the inset. 



[choosing ^12(6) < 0] for b —> tends to K as in 
Reversing the procedure, one thus has found a highly 
nongeneric perturbation splitting an M = 4 JB into two 
M = 2 JBs, which do not undergo further splitting. Note 
that the family Tt(b) is not of the form (|3.1|l : it is im- 
possible to obtain this particular structure with only a 
first-order correction AH, if the latter corresponds to 
AK only. See Sections 6.2 and 6.3 (final paragraph) 
in Ref. 0], where the existence of off-axis JBs in a con- 
tinuum model is studied as an open question. 



V. DISCUSSION 

We have extended the eigenvector expansion developed 
in the previous paper to situations where the eigenvec- 
tors merge and thus are incomplete. The Jordan normal 



9 



-0.98 



-0.99 - 



-1.01 



-1.02 
1.31 



-0.98 







(a) 




® 






® 




- 


® 


- 




® 






® 

® 





-0.99 



f -1.00 



1.32 1.33 1.34 1.35 

Re(to) 



(b) 



@ 







® 



© 



® 



basis is then used. Such a basis for a general matrix 
operator 7i (not self-adjoint) is well known, but here its 
properties have to be considered together with the bilin- 
ear map. The bi-orthogonality of the Jordan normal ba- 
sis, most simply expressed by the metric Q2.8fl . has been 
established: part of it being intrinsic and part of it being 
a conventional choice. As a consequence, 7i.. takes the 
symmetric Jordan- type form Q2.9|) . 

The Jordan basis vectors (like the eigenvectors in the 
non-critical case) immediately solve the dynamics. The 
time evolution of these vectors fj >n is characterized by a 
polynomial prefactor in t, as expected when several ex- 
ponential terms with slightly different frequencies merge. 

Perturbations around critical points are particularly 
interesting, a term eA7i shifting frequencies by a frac- 
tional power of e in equiangular directions. Their study 
allows the small-denominator problem associated with 
(fk, fk) — ¥ at criticality to be handled. 

All these concepts have been illustrated by nontrivial 
examples. By treating criticality, this paper complements 
the previous one an d places the familiar concept of 
critical damping into a general framework. As in |]|, 
the entire formalism can be promoted to the quantum 
domain by turning the state vectors into operators. 
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FIG. 5: Eigenvalues for two M — 2 JBs split by feu i— > feu +e, 
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— 10 -4 . Crosses (circles) denote numerical (perturbative) val- 
ues. Only the right-half oj-plane is shown. The nearly equal 
spacing shows that the shifts are oc e 1 ^ 2 . 
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